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ABSTRACT 

We examine a physical process that leads to the efficient formation of gas gi- 
ant planets around intermediate mass stars. In the gaseous protoplanetary disks 
surrounding rapidly-accreting intermediate-mass stars we show that the midplane 
temperature (heated primarily by turbulent dissipation) can reach >1000 K out 
to 1 AU. Thermal ionization of this hot gas couples the disk to the magnetic 
field, allowing the magneto-rotational instability (MRI) to generate turbulence 
and transport angular momentum. Further from the central star the ionization 
fraction decreases, decoupling the disk from the magnetic field and reducing the 
efficiency of angular momentum transport. As the disk evolves towards a quasi- 
steady state, a local maximum in the surface density and in the midplane pressure 
both develop at the inner edge of the MRI-dead zone, trapping inwardly migrat- 
ing solid bodies. Small particles accumulate and coagulate into planetesimals 
which grow rapidly until they reach isolation mass. In contrast to the situation 
around solar type stars, we show that the isolation mass for cores at this critical 
radius around the more massive stars is large enough to promote the accretion of 
significant amounts of gas prior to disk depletion. Through this process, we an- 
ticipate a prolific production of gas giants at ~ 1 AU around intermediate-mass 
stars. 

Subject headings: planetary systems: formation, protoplanetary disks 
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Introduction 



The discovery of a plethora of extra-solar planets around solar-type main sequence stars 
has established that planet formation must be a common process, not a peculiarity of our 
own solar system. As observational techniques for planetary detection have become more 
sophisticated, the discovery domain has expanded to include host stars with a wide range 
of masses. While on the main sequence, intermediate mass stars (stars with 1.5M Q < M* < 
3M Q ) make poor radial velocity (RV) su rvey candidat e s as th ey have few spectral lin e s whi ch 
also tend to be rotationally broadened (iGrifiin et al.l (120001 ) but see iGalland et al.l (120061 )). 
However, once these stars evolve off the main sequence, their relatively cool and slowly 
rotating outer layers make them more suitable candidates for high precision spectroscopic 
studies. Recent RV surveys targeting evolved intermediate mass stars suggest that they 
differ from solar-type stars as planetary hosts in at least two respects. First, the total 
frequency of giant planets (with periods less than a few years) appears be higher around 
intermediate mass stars. Second, the planets have different statistical properties. Their 
semi-major axis distribution is concentrated at 1-2 AU and there is an apparent lack of 
short-peri od (days to months) p lanets, despite observational selection effects favoring their 
discovery (ILovis fc Mayorl 120071 ). 



In this paper, we propose a common explanation for prolific gas giant formation with 
semi-major axes comparable to 1 AU and for the rarity of close-in planets around intermediate- 
mass stars. As it is unlikely that all planets within 1 AU have been engulfed or had thei r 
orbits disrupted by the current expanded envelope of the host stars (I Johnson et al.l 120071 ). 
we attribute both properties to the formation and early evolutionary processes rather than 
to post- main- sequence evolution. 

We begin by examining the physical properties of circumstellar disks which may affect 
the probability of forming g i ant p lanets. In the core-accretion model of planet formation 
(cf. iBodenheimer fc Pollack! (119861 )). the emergence of Jupiter-like gas giants requires that 
a population of solid cores form within a gaseous protoplanetary disk. These cores grow 
through cohesive collisions with planetesimals, with a growth rate determined by the veloc- 
ity dispersion of the planetesimal swarm. The magnitude of this velocity dispersion is set 
by a balance between excitation by gravitational perturbations and damping by gas drag. 
In the gas-rich environment of typical protostellar disk, gas drag dominates so that field 
planetesimals only attain relatively small equilibrium velocity dispersion. As a result the 
most massive protoplanetar y embryos can acces s only those building blocks within their 
gravitational feeding zones (IKokubo fc Idal Il998h . When these embryos have collected all 
the planetesimals within about five times their Roche radius on either side of their orbits, 
their growth stalls. This maximum embryo mass, a function of planetesimal surface density 
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and distance from the central star, is referred to as the embryo's isolation mass (M; so ). Gas 
giants can only form if the embryos' M iso is suffic iently large for t he cores to begin accreting 
gas prior to the depletion of their nascent disks (llda fc Linll2004l ). 



Although the gravity of lunar-mass embryos is adequate to accrete disk gas with tem- 
perature < 10 3 K, efficient dynamical gas accretion is only possible for cor es with masses 
greater than some critical value (M crit ). In a minimum mass solar nebula (IHayashil 1 19811 ) 
with an interstellar grain size distribution, M crit ~ 10M ffi at a semi-major axis a ~ 5 AU 



(IPol 



ack et al 



1996L although this critical mass decreases both with lowered grain opac- 



i ty (llkoma et al 



2000 



( IBodenheimer fc Pollack 



1986 



Hub ickyi et al.l 120051) and with incre ased density of the ambient gas 



Papaloizou fc Terquemlll999l ). Gas giant formation therefore 



requires that the heavy-elements in the disk can be efficiently assembled into massive cores 
with mass greater than M crit . In order to understand the spatial distribution of the gas giant 
planets we must understand how the building blocks of these cores migrate and are retained 
in gaseous disks. 

In protoplanetary disks solid retention first becomes an issue once grains grow beyond 
a few cm in size. In most regions of protostellar disks, the midplane pressure (P m id) de- 
creases with distance from the central star (r) so that the gas is slightly pressure supported, 
resulting in a sub-Keplerian azimuthal velocity. Grains larger than a few cm are decou- 
pled from the gas and move at Kepleria n speeds. Consequently , grains typically experience 
head winds and undergo orbital decay (jWeidenschillingi 119771 ). However if P mid does not 
monotonically decrease with r then immediately interior to a local pressure maximum the 
gas attains super-Keplerian velocities. This motion introduces a tail wi nd on the decoupled 



grains and causes them to d rift outwards towards local pressure maxima (IBryden et al.ll2000 
Haghighipour fc Bossl 120031 ) . 



Solid retention again becomes an issue once planetesimals grow into earth-mass em- 
bryos and tidal interactions with the gaseous disk be come important. Befo re embryos 
are sufficiently massive to open up gaps in the disks (ILin fc Papaloizou! Il986l ). they can 
exchange angular momentum with the gas via their Lindblad and co-rotation resonances 
( IGoldreich &: Tremaindll980l ). A geometric bias causes an imbalance between the Lindblad 
resonance s which generally leads to a loss of angular momentum and orbital decay for the 
embryos (jWardl Il986l . I1997D . However, embryos will gain angular momentum through their 
co-ro tation resonances if there is a positive P m id gradient (ITanaka et al.l 12002 ; iMasset et al. 
20061 ). Numerical models which take into account thes e physical effects have reproduced the 



observed M p — a distribution around solar type stars (jlda fc Lmll2008l ). 



Several physical processes can lead to local maxima in P m id- Various authors have ex- 
plored the potential accumulation of grains at transient pressure maxima formed by turbu- 
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lent fluctuations ( Johansen et al. 2006 ; Fromang fc Nelson] bood ) or spiral waves ( Rice et al. 



20061 ). These mechanisms, while likely extremely important for forming planetesimals at a 
large range of radii, are still quite "leaky" as a significant fraction of the solid material sim- 
ply undergoes a slightly slower random walk towards the central star. Howeve r, longer lived 
press ure maxima may also exist due to large scale changes in the disk viscosity (IKretke fc Lin 
2007h . 



We expect r adial variations in viscosi ty if turbulence caused by the magneto-rotational 
instability (MRI; iBalbus fc Hawleyl (119911 ) ) is the primary mechanism for transporting an- 
gular momentum. These variations result from changes in the ionization fraction at different 
radii in the disk since free electrons are needed to couple the gas to the magnetic field. The 
disk is thermally ionized in the hot inner regions, but further out stellar x-rays and diffuse 
cosmic rays ionize only the surface lay ers, resulting i n a viscously active turbulent surface 
sandwiching an inactive "dead zone" (jGammid Il996l ). At the critical radius marking the 
inner edge of the dead zone (a cr it), the effective viscosity decreases with increasing distance 
from the central star. In a quasi-steady state situation (expected to develop rapidly in the 
inner regions of the disk) this decrease in viscosity leads to a local increase in the magnitude 
of T, g and hence of P mi d with radius. This disk structure provides a promising barrier to the 
orbital decay of both boulders and embryos. The radial location of a crit depends on the stellar 
mass and the mass accretion rate which we argue explains the observed differences between 
the statistical distributions of planets around solar-type stars and around intermediate mass 
stars. 

In this paper, we present a model for the formation of planets at the inner edge of the 
dead zone and argue why this process is more relevant for intermediate mass than for solar 
mass stars. In §|2]we describe our quantitative nu merical model for the evolution of solids in 
the disk, based upon the work of iGaraudl (120071 ). An important aspect of this model is the 
location of a crit (the inner edge of the dead zone) which we derive in §2.21 as a function of 
stellar mass and mass accretion rate. In $3] we present the model results for a 2 M Q star and 
estimate how this planet formation mechanism scales with stellar mass. In §H]we summarize 
our conclusions. 



2. Model Description 



In order to assess the probability of forming of gas giants at the inner edge of the dead 
zone, we must calculate the expected isolation mass of the cores (M iso ) at this location 
(acrit). From the work of iKokubo &: Ida! (119981 ). it has been established that the embryo's 
isolation mass is sensitive both to the distance from the host star and to the surface density 
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of solids. For the purpose of computing the efficiency of solid retention we calculate the 
dynamical ev olution of the solids in the disk using a modification of the numerical scheme 
developed by lGaraudl (120071 ) (hereafter G07). For full details we refer the readers to G07 but 
for reference the salient points are described here. In this model we assume that particles, as 
a result of a collisional cascade, maintain a power-law size distribution in which the number 
density of particles of size s goes as dn/ds oc s -3 5 , with sizes ranging from a fixed s m i n to 
s max and where s max is allowed to vary with time and distance from the central star. This 
simplifying assumption allows us to completely describe the evolution of the gas and solids 
in the disk by the total gas surface density (S g (r, £)), the solid and vapor surface densities 
of the different species in the disk (S pi (r, t) and S Wji (r, t) respectively), and the maximum 
size of the particles s max (r, t). 

In the following sections we describe how we have modified the G07 scheme. In §2.14 we 
describe how we account for the presence of an evolving dead zone. In particular we describe 
in §2.21 how we calculate the location of the inner edge of the dead zone (a cr it) using a simple 
model of the vertical thermal structure. Finally, in §2.3l we model the drag on the gas caused 
by the particles in the limit of large Z = S p /S g , an effect important near a cr ; t where solids 
are seen to accumulate. 



2.1. Viscosity at a function of r and t 

The original G07 model uses the standard a viscosity prescription (z/(r) = a c ^c s {r)h{r)) 
where c s is the midplane sound-speed, h is the disk scale-height (h = c^id^ 1 ) and a c s is 
constant. In this paper we consider a eff = a eff (r, t) to model the presence of an evolving 
MRI-dead zone. In fully MRI-active regions (eg. interior to a cr i t ) we assume that a c s = «mri, 
while exterior to a cr i t , a e g is modified to take into account the lower viscosity of the dead 
zone. In §2.21 we describe in detail how we calculate the location of a cr it. 



Exterior to a cr i t the disk is ionized by x-rays and cosmic rays, so only the surface layers 
are MRI active. The column density of this active layer (£a) is strongly dependent on the 
ionizing source (ie. the ability of cosmic rays to penetrate the stellar magnetosphere) and 
on the recombination rate, which is dominated by the amount of small grains. Due to these 
uncertainties, and the fact that the exact evolution of the outer parts of the dead zone are 
relatively unimportant in these calculations, we follow previous studies and assume th at 



100 g cm at all radii (outside of the thermally ionized region) (eg. iGammid (119961 )). 



Additionally, we assume that there is some amount of angular momentum transport in 
the dead zone. This transport may be due to the propagation of MRI-driven waves into 
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the laminar dead zone ([Fleming fc Stond (120031 ); iTurner et al.l (120071 )) or to a mechanism 
unrelated to the MRI. This motivates the following prescription in regions beyond a cr i t 



otes(r,t) 



2aMRlSA+aeff,dead(S fl ~2E j4 ) 



E s (r,t) 



«MRI, 



if S g > 2S A , 
otherwise, 



X 



where a e s dead is the effective viscosity in the dead zone. 



2.2. Location of a crit 

As the degree of thermal ionization is a sensitive function of temperature, we must 
calculate the vertical structure of the disk in order to find the location of a cr it, the inner edge 
of the dead zone. In the innermost regions of the disk viscous heating dominates the energy 



budget, and stellar irradiation can be neglected (IGaraud fc Linl 120071 ). 



We assume that the disk is in hydrostatic equilibrium 

^ = - Pg n 2 K z, (2) 

and is viscously heated so that the vertical energy flux (F) is described by 

Here we consider that z/(r, z) is a function of both height and radius in the disk and assume 
that the orbital frequency can be approximated by the Keplerian frequency. 

We also assume radiative energy transport so that 

_ 4ac T 3 dT 
3 Kp g dz ' 



(4) 



where k is the Rosseland mean opacity. For the temperature range of interest (400-1, 200K) 
the dominant opacity sources are silicate and iron grains. We adopt a grey opacity approx- 



i matio n in which n = 1 cm g . This approximation is consistent with the iFerguson et al. 



( 120051 ) dust opacities. We solve only the optically thick region of the disk and use photo- 
spheric boundary conditions at z = z e , namely P(z e ) = (2/3)Q 2 K z e / k and F(z e ) = aT(z e ) 4 . 
Finally we assume symmetry about the midplane. 

In order to parametrize the variation of the viscosity with height above the midplane we 
follow the results of 3D MHD simulations which demonstrate that, in fully developed MRI- 
turbulence, the shear stress(w = up g r(dfl/dr) = (3/2)z/p 9 f2 K ) is approximately constant 
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with height (jMiller fc Stondl2000l ). Based on these numerical results, we model the viscosity 
at a given radius as 



2 qPjs 

3 p 9 (z)n K 
0, 



otherwise, 



(5) 



where z v ~ 2/i. Note that this differs f rom the more co mmon 2D parametrization sometime 

{z)VL^ . Viscosity 



(e.g. 



CannizzJ (jl992l )) where v a _ 



ac: 



referred to as the "aP-formalism' 
described by equation ([5]) increases with height, which leads to a cooler mi dplane more 



consi stent with multi-dimensional MHD simulations than the aP-formalism (IHirose et al. 
2006h . 



Equation (jSJ) is applicable to MRI-active regions where the partially ionized gas is well 
coupled to the turbulent magnetic field through the entire thickness of the disk. As long 
as the midplane temperature is high enough to sufficiently ionize the gas this structure 
is self-consistent. We can therefore use it to determine a cr it, the outermost radius where 
this condition is satisfied. In practice we combine the equations for the vertical structure 
with the Saha equation to calculate the location at which the midplane just satisfies the 
ionization criteria of x e > 10~ 12 (where x e is the fractional numbe r density of electr ons) 
which corresponds to T ~ 1000K for typical disk midplane densities (lUmebayashil Il983l ) . 



In Figure [2] the curves show the location of the inner edge of the dead zone as a function 
of stellar mass and accretion rate assuming that in the active region a = 10~ 2 . These 
theoretical curves can be approximated by 



0.77 



M 



4/9 



X 



1.4 x lO" 7 M yr- ] 
a \ - 1 / 5 ( k d 



2M Q 



1/3 



10 



-2 



1 cm 2 g 



1/4 



AU, 



where the quasi-steady-state approximation for the mass accretion rate 

POO 

M = 3tt / pgudz = 47raP m id^K 1 -2 l , 



(6) 



(7) 



was used to eliminate P m \<\- We find that this value for a cr i t varies significantly from the 
~ 0.1 AU often quoted (ie. iGammid (119961 )) if the stellar mass or the mass accretion rate 
are large. By using the relationship 

M = 37tz/S 9 = 37ra cff c 2 ft K % (8) 



( jPringld Il98ll ) we can relate the a in equation (jSJ) to the the vertically averaged a e s when 
the disk is fully MRI active and find that omri ~ ol. 
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Observationally, M appears to be correlated with M#, although the exact relationship 
is both uncertain and shows significant scatter. For reference the symbols in Figure 2 in- 
dicate measurements of mass accretion rates onto you ng stars. Th e solid points are from 
observations of the young cluster (< 1 Myr) p-Oph by iNatta et al.l (120061 ) and the dashed 
line shows the best fit to their data which is 



M ~ 4 x lO" 8 (M,/M ) L8 M yr- 



(9) 



As the mass accretion rate for the higher mass stars in the Natta et al. sample is dominated 
by a single object, we have also plotted (as open points) similar data from older, hetero- 
geneously distributed inte rmediate mass stars with estimated ages ranging from 1-10 Myr 
(IGarcia Lopez et al.l 120061 ) . As suggested in this study, these systematically older stars may 
have had higher mass accretion rates consistent with the extrapolation from p-Oph, when 
they were younger. 



Using the best fit relationship from p-Oph implies that a cr i t oc M} A . IClarke fc Pringld 



(120061 ) suggest that the correlation between stellar mass and mass accretion rate may not 
be this steep due to observational biases. If instead we use their estimation that M oc M*, 
then a crit oc M^ 8 , which is a slightly less sensitive but nevertheless increasing function of 
M*. Therefore, in the light of these the uncertainties, we will simply assume that a cr i t oc M*. 



2.3. Gas-Particle Feedback 



In a thin disk, the global evolution of the gas surface density £ 9 is determined by the 
equation 

r)Y,. 1 F> 

0. (10) 



<9X„ 1 c> 
at r ar 



In a standard viscous accretion disk the radial velocity (u r ) of gas equals u v where 

3 d 



u v = 



r 1 / 2 E 9 dr 



(11) 



In our analysis, we consider the possibility that the radial velocity of the gas u r is not only 
determined by the viscosity but also by momentum transferred via drag between the gas 
and the solid partic les. We calculate the eq uation of motion for a parcel of gas and dust in 
a similar fashion to iNakagawa et al.l (119861 ). However, instead of assuming a single particle 
size for the solids we consider the power-law distribution (dn/ds oc s~ 3 ' 5 ). The equation of 
motion for a given particle of size s is 



dV(s) 
dt 



tAs 



'V(s) 



rn GM * 



(12) 
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where r s (s) = sp s /(p g c s ) is the stopping time for a particle of size s and density p s (in the 
Epstein regime, see G07 for a full description), where V(s) is the size dependent particle 
velocity and U is the gas velocity. 

Assuming that deviations from Keplerian orbital velocity are small, we write V(s) = 
v r (s)r + (rfix + v<j,(s))4> and U = u r f + (rfix + U(j>)4")- The linearization of equation (TT2l) 
yields 

dv I 

-(v r (s) - u r ) + 2Cl K v^(s), 



Of 
dt 



T S {S / 

■y—r (v^s) ~ Uf) ~ -Q K V r (s), 



(13) 



The particles and gas adjust to a steady motion with respect to each other within a few 
stopping times, so we solve for the steady-state solutions only. Multiplying by the particle 
mass and integrating over the whole size distribution function yields 



mm 

max 



dn 1 ^ . .dn . . 

-m(s)— —{v r {s) — u r ) + zil-^misj—vAs) 

ds 



ds t s (s 
dn 1 



I dfi 

-^IsJ — — — (^(s) - tty) - -Q K m{s)—v r {s) 



ds 



ds 



0. 



0. 



(14) 



The momentum lost by the solids is acquired by the gas, so that the steady state gas-dynamic 
can be described by 



m(s 



mm 
5 max 



dn 1 

ds tJs) 



[v r (s) - u r )ds + 2£l K p g (u ( j } + nv K ) 



m(s)^— j— (^(s) - u^ds - ~Q K p g (u r - u v ) 
as T s {s) l 

where n is the non-dimensional pressure gradient 

lh 2 d\nP 

V = 



0. 
0. 



(15) 



(16) 



2 r 2 Q\ nr 

and vk = rQ K . In most regions of the disk n is positive due to the negative pressure gradient, 
but its sign will change around pressure maxima. In the absence of particles, the gas will 
have an azimuthal velocity (relative to the Keplerian motion) of = —tjvk and a radial 
velocity of u T = u v . 



Combining these equations and defining the mass-weighted average particle velocities 



as 



v<t> = — 
1 



mm 
max 



m(s)-^v <t> (s)ds, 

m(s)—v r (s)ds 
ds 



(17) 



-10- 



allows us to write 



Pa 

Combining this with equation (fT51) we can solve for 

u r + 2Q K r s (s)u l p 



P P - 

U<j> = — -v^, 
Pg 

u r = u v — —v r . (18) 



v r (s) 



l + 0|r s ( S ; 
- (l/2)Q K r s (s)u r 



= 'uw • (19) 

Using the definition for the Stokes number St(s) = fiK r s(s)/(27r) and taking the mass- 
weighted integral over all particle sizes yields 

v r = I (V27rSt max ) Mr + 2J (V2vrSt max ) 

v<j> = I (V 27rSt max) - ^ J (V 27rSt max) u r , (20) 

where I and J are the same as equation (52) in G07 (reproduced in Appendix A for reference) 
and St max = St(s max ). Combining equations ([TBI and f[2"U|) yields the final steady-state 
velocities for the particles and the gas where 1 = 1 (y/2nSt max ) , J = J (v / 27rSt max ) , and 

X = Pv/Pg- 



Ur 



[I + X (I 2 + J 2 )]u„ - 2J V v K 
1 + 2 X I + X 2 (I 2 + J 2 ) ' 

1 Uu u + 2[l + x (i 2 + J 2 )}yv K 

2 I 1 + 2 X I + X 2 (I 2 + J 2 ) 
(1 + xl)u v + 2Jxvvk 



(21) 



1 + 2 X I + X 2 (I 2 + J 2 ) 
1 

2 



and 



- 2(1 + xI)vvk 
l + 2 X I + xV 2 + J 2 ) 



(22) 



In the gas dominated limit (x — > 0) these equations reduce to those in G07. If St max 
is large and St max ^> x 2 ■, the solids will be decoupled from the gas and will not migrate 
significantly (v r = 0) while the gas evolves viscously (u r = u u ). In the limit St max <C 1, 
the grains are well-coupled to the gas so there will be little relative motion between the two 
(u r = v r = u u /(l + x))- 
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3. Model Results 

With these modifications to the G07 prescription we calculate the evolution of a disk 
around a 2M star including an MRI-dead zone (qjmrj = 10~ 2 , a e fr,dead = 10~ 3 ). The initial 
disk has a surface density profile of 

E » = E »^nu)" exp (^ (23) 



'0 



where Ro = 30AU, and So = 10 4 g cm -2 has been chosen such that the quasi-steady state 
accretion rate M from equation ([8]) is 1.4 x 10~ 7 M yr^ 1 when a e g = a e ff,dead- For the chosen 
value of a e ff 5 dead the total disk mass is 0.1M*. This relatively massive disk is still stable 
according to the Toomre criterion at all radii. For simplicity we only track one species 
of solids, a generic refractory material which we take to be a combination of silicates and 
metals, materi als which are assumed to sublimate at 1500K. Using the solar composition from 



Loddersl (120031 ) we begin with a dust-to-gas ratio of 0.005 for these refractory materials. We 
do not track volatile ices as they will not contribute to the solid cores formed in the hot 
regions of interest. For other model parameters we use the values presented in the fiducial 
model of G07. 



3.1. Accumulation of Solids 

We first study the accumulation of solids near the edge of the dead zone and emphasize 
the importance of feedback in terms of momentum exchange between the solids and the gas. 
We present two runs, the first neglecting momentum transfer from the dust to the gas and 
then including this feedback. 

The top and middle panels of figure [3] show the evolution of the gas and solids neglecting 
the feedback of the dust on the gas. For reference, figure H] shows a e ff at 10 4 years and at 10 5 
years to indicate the location of the dead zone. Early on the dead zone extends from 0.9 to 
20 AU, but as the disk evolves the dead zone shrinks. Within the first 10 4 years the gaseous 
disk adjusts to a quasi-steady-state profile in the inner regions. The surface density of the 
gas in the fully MRI-active inner region (r < 1 AU) is reduced by a factor of a e ff,dead/aMRi 
compared with that in the dead zone (r e [1, 10] AU), thus creating a local pressure maximum 
near 1 AU. In this calculation neglecting feedback, the solids accumulate in a very narrow 
ring corresponding to this pressure maximum. After only 5.3 x 10 5 years the largest-size 
body at this location reaches 5M ffi and the core thus formed is expected to continue to grow 
significantly due to accretion of gas, not included in this numerical calculation. Crucially, in 
this model tidal interactions will not affect the orbital evolution of the core as the effect of 
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Type I migration would be to keep the core at the pressure maximum (IMasset et al.ll2006l ). 
It is also interesting to note that when feedback is neglected, virtually all solids are trapped 
at the pressure maximum and negligible amounts of heavy elements accrete onto the star. 
This effect is associated with the clear spatial separation of the pressure maximum and the 
sublimation line (here at ~ 0.3 AU) and would appear to predict that intermediate-mass 
stars should be strongly depleted in refractory elements. However we now show that this is 
in fact unrealistic as the large build up of material makes it necessary to include the transfer 
of momentum from the solids to the gas. 

The bottom panels of figure [3] show the evolution of the same initial disk including this 
momentum feedback. The early evolution is essentially similar for t < 10 4 yrs. However, 
once the amount of solids in the inner edge of the dead zone has increased by an order of 
magnitude, it begins affecting the gas properties (compare the middle and bottom panels of 
fig[3]). The mass accretion rate onto the star decreases as the gas receives angular momentum 
from the solids. The decrease in M reduces the surface density and thus the midplane 
temperature, causing a cr it to move inwards. As a cr i t moves inwards the particles of different 
sizes respond to this shift in location at different rates, smoothing the sharp peak in the solid 
surface density and allowing the trapped gasses to accrete onto the central star. In this way 
both M and the position of a crit oscillate with time as seen in figure [5j It is also interesting 
to note that these oscillations allow heavy elements to accrete onto the star since when a cr i t 
moves outwards it leaves some solids interior to the "trap" of the pressure maximum. Thus 
intermediate-mass stars will not, in fact, be significantly depleted in refractory elements. In 
these simulations we find that the total amount of solids trapped near the pressure maximum 
is of the order of 100 M @ and that the dust-to-gas ratio in this region is of order unity. The 
core growth is slower than in the simulation without feedback as material is deposited over 
a wider range of radii. In 5.3 x 10 5 yrs the largest body has grown to 0.5M e but the amount 
of material in its vicinity suggests that it will continue to grow larger still. 



3.2. Core Formation 



As the core mass increases much beyond 0.5 M® many effects not included in the 
numerical algorithm could affect its growth. For example, it may be important to include 
the effects of tidal interactions since the oscillations of the position of the pressure maximum 
could affect the orbital evolution of the core. Additionally, gravitational interactions between 
solids may start to impact their radial distribution. The cores may scatter each other out 
of the region of interest and radially migrati ng solids will be trapped into r esonances with 
the existing core, stalling radial migration (jWeidenschilling fc David Il985l ). As a result, 
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we do not continue integrating the simulations beyond this point but instead estimate the 
achievable core mass with simpler but robust scalings based on our most important results 
from the previous section: (1) the presence of a large reservoir of solids (> 100M©) around 
a crit and (2) the saturation of the dust-to-gas surface density ratio to Z ~ 1 near a crit (as 
seen in fig. [3]), which is a simple consequence of the saturation of the nonlinear momentum 
feedback between the solids and the gas. 

U sing these ideas , we estimate the maximum achievable core size from the isolation 



mass (jlda fc Linll2005l ) at a crit is 

M iso = 2ir 1 E p a crit br H ( 

^crit 

), (24) 



where r H (a clit ) = a crit (2M iso / (3M*)) 1 / 3 is the Hill's radius and b ~ 10 (lLissauerlll987l ). We 
assume that S 5 is related to M by the quasi-steady-state approximation ([8]), and that E p is 
related to the gas density via Z = S p /S g where Z ~ 1. We find that 



M sn ~ 60 Z 3/2 



11/6 

M 



x ( — ) {YrjK) Cio=5-J (25) 



1.4 x lO- 7 M yr 

1/2 / N -3/2 

M 

Note that this estimate should be considered as a maximum achievable size and that, more 
realistically, cores may only reach a fraction of M iso . Therefore in the following section we 
select the core mass M c = fM lso . The isolation mass of the core at a crit goes as M iso oc 
M n / 6 M* 1/2 . If we assume that M is related to M* as discussed in JTlthen M iso oc Ml to 
Mf. Therefore, super-earth cores form preferentially at the inner boundary of the dead zone 
around more massive stars. 



3.3. Gas Giant Formation 



As our numerical algorithm does not include gas accretion onto the cores we can use 
the estimate of core mass derived in the previous sect ion to look at the potential for further 
growth using the approximations of llda &: Linl (120051 ). The growth of a planet of mass M p 
due to the gas capture can be approximated as 



dt 



T 



where 



T 



To 



Mp 



-3 



(26) 



(27) 
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with To ~ 10 10 years (see llda fc Linl (120051 ). although tq is likely to be smaller in the inner 
regions of the disk where p g is relatively large). As long at the final mass of the planet is 
much greater than the original core mass, the timescale for giant planet formation is then 



t 



giant 



3 



Mc 



(21 



This inference effectively means that in order for a gas giant to form within the lifetime 
of the evolving disk, the core mass must at least be of the order of 10 M e . Using the 
relationship derived in equation (!25l) implies that t g i ant oc M~ 11/2 M*~ 3/2 . Folding in the 
relationship between M and M* yields t g i a nt oc M~ 7 to M~ 12 . This very steep function 
demonstrates that giant planets will not have time to form at the inner edge of the dead 
zone around less massive stars, but will form ubiquitously around higher mass stars with 
higher accretion rates. 



3.4. Asymptotic Mass and Multiple Planet Systems 



The observationally inferred mass of the known planets around stars with M* > 2M & 
is in the range of 2 — 20Mj which is near the upper end of the M p distribution of known 
planets around solar- type stars. In the core-accretion s cenario, the asymptotic mass of gas 
giants is determined by a thermal truncation condition (ILin fc Papaloizoulll993l ) such that 



M p ~ f a {h/r) 6 M* 



(29) 



where the co nstant f a is of the order o f 1 to 10 and depe nds on the detailed thermal structure 
of the disk (IDobbs-Dixon et al.l 120071 ; llda &: Lirj 120051 ). We use equation (Q to determine 
the value of h/r = c s /vk at a crit . The temperature and sound speed at the inner edge of 
partially dead zone is, by construction, independent of M* and M. We find from equations 
(J6]) and (J9j) that, at a cr it, h/r ~ 0.05 during the main phase of disk evolution. From this 
result, we estimate that 

h/r\ 3 ( M, 



0.3/ a 



0.05 



2M 



Mr. 



(30) 



Figure 2 shows that the dependence of M on M* in equation ([9]) is an average relation 
with considerable dispersion. Furthermore, M declines with the stellar age £*. Within our 
model planets with larger M p on longer period orbits can form during the early epochs of disk 
evolution when M is much larger and a cr it is correspondingly further out. However, the rapid- 
accretion phase may be too brief for embryos to reach their isolation mass. Nevertheless, 
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cores with sufficient mass to initiate efficient gas accretion are likely to emerge when M is 
reduced to that approximated by equation (Q and when the timescale for core formation at 
a crit becomes comparable to the disk evolution timescale. Additionally, due to this spread in 
M around stars with similar M* we anticipate there will be a wide distribution of planetary 
masses, albeit the mean value of M p should vary with M*. 

In disks with protracted high-M evolution, the first generation of gas giants forms 
rapidly. Once the first planet has formed and grown large enough to open a gap, the outer 
edge of this gap provides another pressure maximum ca pable of trapping solids and promoting 
the formation of the next planet (IBryden et al.l (120001 ) ). We anticipate a prolific production 
of multiple planetary systems for intermediate-mass stars. 



4. Summary and Discussion 

While planets are likely to form by the same basic mechanisms regardless of the envi- 
ronment in which they form, the properties of their host stars and the detailed structure 
of their nascent disks will strongly affect the statistical outcome of the formation process. 
Observations hint that planets may form systematically more efficiently around intermediate 
mass stars, and that there is a statistically significant lack of giant planets on orbits with 
semi-major axes much smaller than 1 AU. In this paper we propose a mechanism which 
forms giant planets preferentially around intermediate mass stars with radial distributions 
roughly consistent with these observations. In this model the gaseous protoplanetary disk 
evolves due to MRI-driven turbulence, creating a pressure maximum at the inner edge of the 
dead zone (a crit ) which traps solid material. In order for the cores formed at this location to 
grow large enough to seed giant planets, the inner edge of the dead zone must be sufficiently 
far from the host star. We demonstrate that, as a cr i t is roughly proportional to M*, this 
condition is only likely to be met around intermediate mass stars. 

The amount of solids which accumulates near a crit is limited by momentum feedback 
on the gas by the solids. This is interesting as it means that the surface density of solids 
at a crit depends on the gas surface density and not on the initial fraction of solids, except 
in extremely metal poor disks. This suggests that while there may be a critical metallicity 
required in order to form planets by the mechanism described in this paper, beyond that 
critical value the frequency of giant planets should depend only weakly on stellar metallicity. 

The prolific production of gas giants near a cr i t can also promote the emergence of addi- 
tional gas giants at larger distances from the same host stars. We may expect the fraction 
of intermediate-mass stars with multiple Jupiter-mass planets is likely to be larger than that 
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around solar-type stars. Nevertheless, we anticipate the peak in the planets' semi-major axis 
distribution to be around 1 AU. This corresponds to the location of the original pressure 
maximum at the inner edge of the dead zone. Quantitative verification of this expectation 
requires population synthesis which will be carried out and presented elsewhere. Observa- 
tional confirmation of this peaked period distribution will provide clues and constraints on 
the outstanding issue of magnetic turbulent transport in protostellar disks. 
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A. Functions I and J 

The functions I and J used in equation (120]) and (122]) are 

= ^[/i(x) + / 2 (x)], 

n 1 ( x 2 + xV2 + l 
[X) 2 11 y x 2 -xy/2 + l 

fz{ x ) — arctan + 1 J + arctan (xy/2-lY (Al) 

In the limit of x -C 1, I — 1 and J = x 2 /3, and in the limit of x 1, / = J = V2n/ (Ax) ~ 
1/x. 
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Fig. 1. — A comparison of the disk structure resulting from the viscosity prescription from 
eq. (JSJ) (solid curves) to disks with the traditional "a-P formalism" (dashed curves). The 
panel (a) shows the accretion stress, panel (b) the viscosity, and panel (c) the temperature. 
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Fig. 2. — The curves show the position of the inner edge of the dead zone as a function of 
stellar mass and mass accretion rate. Symbols represent observations of mass accretion rates 
for stars of various masses (see t ext fo r details). The dashed line shows the best-fit for the 
p-Oph cluster from iNatta et al.l (120061 ). 
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Fig. 3. — Evolution of the surface density of solids (solid curve), gas (dashed curve) and 
heavy element in vapor form (dotted curve). The middle panels show the evolution neglecting 
feedback of the dust on the gas while the bottom panels include feedback. 
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Fig. 4. — The calculated a cS as a function of radius at 10 4 and 10 5 years. The inner edge of 
the dead zone, a crit is indicated by the dotted line. 
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Fig. 5. — The mass accretion rate onto the star as a function of time. The solid curve 
indicates the mass accretion rate when feedback is included while the dashed curve neglects 
feedback. 



